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INTRODUCTION 


In  a  number  of  recent  papers,  in  particular  by  Martin,  Symonds  and  co¬ 
workers  [1966-1969]  approximate  solutions  have  been  sought  to  the  response  of 
plastic  structures  subjected  to  dynamic  loading,  and  general  theorems  for  the 
rational  choice  and  the  evaluation  of  the  (average)  reliability  of  such  approx¬ 
imations  have  been  established. 

However,  the  applications  of  these  theorems  have  been  mostly  limited  to 
impulsive  loading  problems  (i.e.  problems  in  which  an  initial  set  of  velocities 
is  imposed  on  the  structure ) .  When  other  loads  were  allowed  to  act  on  the  struc¬ 
ture,  they  had  to  be  below  the  static  collapse  value  [Oien-Martin ,  1965]. 

Only  very  recently  the  procedure  has  been  extended  to  single-pulse  loading 
problems,  i.e.  problems  in  which  a  high  load  is  suddenly  imposed  on  the  struc¬ 
ture,  and  then  either  is  suddenly  removed  after  a  given  time-interval  or  decays 
more  or  less  slowly,  tending  to  zero  as  time  increases.  (Such  sort  of  loading 
can  be  used  for  instance  to  approximate  blast  loading . )  It  has  been  shown  that 
a  rational  approximation,  although  less  satisfactory  than  for  the  impulsive  load¬ 
ing,  is  possible  in  this  case  also,  at  least  within  certain  limiting  assumptions 
that  may  be  considered  satisfied  in  most  actual  instances  [Augusti,  Martin, 
O'Keeffe,  1969], 

Independently, very  similar  results  were  obtained  by  Kaliszky  [1968]*  who,  how¬ 
ever,  lacked  a  general  measure  of  the  approximation,  although  in  the  several  ex¬ 
amples  in  which  he  compared  his  approximate  results  with  the  exact  ones  a  very 
good  agreement  was  shown. 


The  writer  is  grateful  to  Dr.  Kaliszky  for  letting  him  have  an  interim  copy  of 
the  full  text  of  this  paper. 
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As  an  essential  feature  of  all  the  works  quoted  above,  the  loading  program 
is  such  that  the  motion  of  the  rigid-plastic  structure  lasts  only  a  compara¬ 
tively  short  time,  after  which  the  structure  comes  to  rest  in  a  permanently  de¬ 
formed  configuration.  Essentially,  no  strain  reversals  take  place  during  the 
response . 

In  another  important  class  of  dynamic  problems,  a  structure  may  be  subjec¬ 
ted  to  a  comparatively  large  number  of  repeated  and/or  reversed  loads,  as  for 
instance  during  most  strong  earthquakes.  Problems  of  this  type  have  been  so 
far  treated  almost  exclusively  by  numerical  means ,  especially  in  the  plastic 
range  of  structural  behavior,  and  no  general  theorems  have  been  established. 

The  study  of  which  this  report  presents  the  first  results, was  undertaken 
with  a  view  of  eventually  filling  this  gap  and  furnishing  procedures  for  the 
dynamic  analysis  of  plastic  structure  under  any  loading  history.  Therefore, 
rather  than  solving  particular  problems,  the  attention  is  focused  on  recogniz¬ 
ing  general  characteristics,  outlining  avenues  of  approach,  and  so  on. 

Unfortunately,  it  has  been  proved  impossible  to  apply  to  general  loading 
histories  the  simple  approximate  solutions  suggested  for  impulsive  and  single¬ 
pulse  loads,  while  keeping  a  rational  measure  of  the  approximation  involved. 

A  substantial  portion  of  the  report  deals  in  particular  with  periodic  load¬ 
ing  of  structures  which  behave  symmetrically  with  respect  to  the  sign  of  the 
loads;  it  will  be  possible  to  derive  some  rather  general  results  on  the  rigid- 
plastic  response  in  this  case ,  without  any  qualifications  of  the  structural 
type.  However,  the  symmetry  assumptions  are  particularly  realistic  in  the  case 
of  rectangular  frames  subjected  to  horizontal  loads,  and  specific  results  per¬ 
taining  to  these  structures  are  also  presented,  including  multi-degree  of  free¬ 
dom  response . 
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The  whole  study  has  so  far  been  limited  to  rigid-perfect ly  plastic  be¬ 
havior  only:  some  aspects  of  this  assumption  are  briefly  examined  in  the  Con¬ 
clusions  , 

1.  SINGLE-MECHANISM  RESPONSE  OF  A  RIGID-PLASTIC  STRUCTURE  TO  DYNAMIC  LOADING 

Consider  a  rigid-plastic  structure,  initially  at  rest,  and  subjected  to 
variable  loads  ,  proportional  to  one  time -dependent  parameter  £  only: 

P4(t)  =  C(t)P£  (1.01) 

where  P£  is  a  (fixed)  reference  value  of  the  A**1  load  (force  or  couple,  con¬ 
centrated  or  distributed).  For  convenience  of  the  following  discussion,  it 
will  be  assumed  that  P£  has  such  dimensions  that  ?(t)  is  an  acceleration. 

Within  this  definition,  the  "loading"  can  be  constituted  by  movements  of 
the  external  restraints  of  the  structure ,  in  such  a  way  that  the  distribution 
of  their  accelerations  is  constant  with  time.  In  this  case,  5(t)  measures 
directly  the  acceleration. 

Just  to  fix  ideas,  consider  the  very  simple  example  in  Fig.  1.1:  it  is 
evident  that  an  applied  horizontal  force  £P  and  a  horizontal  ground  acceler¬ 
ation  Ug  =  -£  are  equivalent  if  the  reference  value  is  chosen  such  that  P=m 

This  section  will  deal  with  the  case  that  the  mechanism  (i.e.  the  shape  or 
mode)  of  the  plastic  response  of  the  structure  is  uniquely  determined.  There¬ 
fore  ,  the  fields  of  displacements  and  strains ,  of  velocities  and  strain  rates , 
of  accelerations  and  strain  accelerations,  are  also  uniquely  determined  through 
out  the  structure  to  within  a  single  multiplying  factor,  which  depends  on  time 
and  can  be  taken  respectively  equal  to  the  scalar  displacement  u  ,  velocity  u 
and  acceleration  u  of  an  arbitrary  point  of  the  structure . 
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Let  us  indicate  by  k+  and  -k  the  values  of  £  (respectively,  posi¬ 
tive  and  negative),  corresponding  to  static  collapse  of  the  structure.  It  is 
immediate  to  recognize  that  the  equations  of  (rigid-plastic)  motion  of  the 
structure,  with  the  given  definitions  and  assumptions,  reduce  to 


u  = 

£  - 

kT 

when 

u  *  0 

u  = 

£  + 

k“ 

when 

u  $  0 

(1.02) 

u  = 

0 

when 

• 

u  = 

0  and 

-k~  $  £  S  k+ 

For  instance,  the  motion  subsequent  to  an  impulsive  loading  in  the  positive- 
£  direction*  shows  a  constant  deceleration  and  a  linear  velocity 

ii  =  -k+  ,  u  =  uq  -  k+t  £  0  (1,03) 

and  a  total  response  time 

t  =  (1.04) 

f  k+ 

where  uq  is  the  (imposed)  initial  velocity  at  t  =  0  .  At  all  times  t  £  t^.  , 
the  third  Eq.  (1.02)  is  verified,  and  the  structure  remains  at  rest.  If  the 
response  mechanism  varies  in  time,  the  above  considerations  apply  to  the  "mode 
approximation11  [Martin  and  Symonds ,  1966]. 

1 . 1  Periodic  Loading;  Case  I 

The  assumption  of  symmetric  structural  behavior  with  respect  to  positive 
and  negative  loads  implies 

k+  =  k"  =  k  (1.05) 

* 

In  a  £  vs.  t  diagram,  the  impulsive  load  can  be  represented  as  a  delta 
function. 
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It  will  now  be  further  assumed  that  the  loading  is  periodic-symmetric,  i.e. 


£(t  +  T)  =  £(t) 
5(t  +  ~)  =  -  £(t) 


(1.06) 


where  T  is  the  period .  An  example  of  loading  vs.  time  diagram  of  this  type 
is  shown  in  Figs.  1.2a  and  1.3a,  where  the  accelerations  given  by  Eqs.  (1.02) 
(for  structure  initially  at  rest)  are  shaded. 

In  fact,  the  structure  does  not  yield  unless  the  loading  function  crosses 
either  of  the  C  =  ±k  lines .  Without  lack  of  generality  we  can  assume  that 
S(t)  crosses  the  £  =  +k  line  first,  say  at  time  t  =  t^  ,  and  is  larger 
than  k  for  a  finite  time  interval  thereafter.  According  to  the  previous  dis¬ 
cussion,  the  acceleration  u  is  given  by  the  first  Eq.  (1.02)  as  long  as 

ft 


u  = 


u  dt  i  0 


(1.07) 


't 


11 


i.e.  for  times  t  comprised  in  the  interval 


"12 


(1.08) 


where  t^.  has  been  already  defined,  and  t^  is  determined  by 


u(ti2) 


'12 


'll 


ft 


u  dt  = 


12 


'll 


(£  -  k)  dt  =  0 


(1.09) 


For  the  simple  load-time  diagram  sketched  in  Figs.  2a  and  3a,  Eq.  (1.09)  implies 
the  equality  between  the  two  shaded  areas  A^  and  A^  • 

At  time  t  =  t^  »  a  discontinuity  of  u  takes  place .  Either  of  two  con¬ 
ditions  can  be  realized,  as  exemplified  in  Figs .  1.2  and  1.3  respectively:  namely, 
either  ^(t.-)  *  -k  (Fig.  1.2)  or  £(t12)  <  -k  (Fig. 1.3).  In  the  first  case, 
which  will  be  indicated  as  Case  I  from  now  on , 
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u(t )  =  u(t )  =  0  for  t12  <  t  i  t13 


At  t  =  t^3  ,  determined  by 

e<*i3>  ■  -k 


or,  equivalently,  by 


'13 


+ 


1 

2 


(1.10) 


(1.11) 


(1.11’) 


a  phase  of  motion  with  negative  velocity  starts,  and  so  on. 

A  periodic  motion  with  period  T  is  thus  activated,  as  qualitatively 
illustrated  by  Fig.  1.2: 


12 

"  "I! 

=  *14 

-  ha 

=  *n2 

=  *114 

“  ^3 

=  A 

=  A 

=  A 

=  A 

=  A  ^ 

=  A  ^ 

=  A 

11 

12 

13 

14 

nl 

n2 

n3 

n4 

T 

2 


=  max|u| 


(n  =  1,2,3. . . ) 


(1.12) 


The  displacement  u  (periodic  but  not  symmetric)  varies  between  0  and 
a  positive  value  Au  ,  equal  to  the  area  of  each  of  the  u-pulses: 

rt  _  ft 


0  *  u(t)  ?  Au  = 


n2 


nl 


u  dt  =  - 


n4 


n3 


u  dt 


(1.13) 


It  is  of  interest  to  note  that  the  total  dissipated  plastic  work  (which 
could  cause  failure  of  the  structure  because  of  large-strain  fatigue)  is  pro¬ 
portional  to  the  accumulated  absolute  values  of  the  displacements*.  Then  , 


By  assumption,  the  mechanism  is  unique;  therefore  the  strain-vectors  do  not 
change  direction  during  the  motion.  When  they  change  sign,  the  associated 
stress  points  jump  to  points  symmetrical  with  respect  to  the  origin  of  the 
stress  axes,  which  by  assumption  is  a  center  of  symmetry  for  the  yield  sur¬ 
face  .  It  will  be  seen  in  Section  3  that  this  behavior  implies  that  each 
stress  point  is  on  a  principal  direction  of  the  relevant  yield  surface. 
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after  n  cycles,  the  dissipated  plastic  work  is  given,  to  within  a  positive 
factor,  by 

u  =  2nAu  (1.14) 


Therefore ,  Au  measures  both  the  maximum  absolute  displacement  and  the 
plastic  work  dissipated  in  each  cycle. 

1 . 2  Periodic  Loading;  Case  II  and  "Steady  State"  Response 
On  the  contrary,  if 


C(ti2)  <  -k 


(1.15) 


the  case  of  Fig.  1.3  is  verified,  which  will  be  indicated  as  Case  II  in  the  fol¬ 


lowing.  The  initial  instant  of  the  first  negative-velocity  phase,  t^  ,  coin¬ 


cides  with  the  final  instant  of  the  first  positive-velocity  phase ,  t 


12  ' 


t  =  t  >  t  +  — 
12  13  11  2 


(1.15') 


the  acceleration  is  given  by  the  second  Eq.  (1.03)  for 


t  <  t  <  t  =  t 
13  14  X21 


(1.16) 


where  t.^  is  defined  by 


u(tm>  * 


’14 


'12 


u  dt  = 


14 


12 


(£  +  k)  dt  =  0 


(1.17) 


(i.e.,  in  the  simple  example  in  Fig.  1.3a,  by  the  equality  of  the  areas  A.^  and 
A  )  and  t^^  is  the  starting  time  for  the  second  cycle  of  motion. 

It  can  be  proved  that ,  provided  k  >  0  , 


*21  '  'll  »  T  •  i-*’  >  k  ; 


(1.18) 


-8- 


therefore,  beyond  ,  the  acceleration  is  again  given  by  the  first  Eq.  (1.02) 

and  the  motion  continues  as  sketched  in  Fig.  1.3. 

The  following  inequalities  can  also  be  easily  proved  (the  relevant  quanti- 


1. 

3a,  b,  c,  all  as  positive 

,  > 

A  >  .... 

>  .... 

11 

21  ' '  nl 

< 

<  .... 

13 

ft23  *  ’  *  ‘  n3 

> 

21 

A.  „  ;  . . . . ;  A  ,  > 
13  ’  nl 

A(n-1)3 

(1.19) 


Au..  >  fiu.,  >  . . . .  >  Au  >  .... 
11  21  nl 


4u12  *  4u22  ‘  ••••  '  4un2  <  " 


4u21  ”  4u12  ;  ••••  ’  4unl  =■  4u(n-l)2 


(1.20) 


It  is  then  intuitive  that  the  motion  u(t)  tends  to  a  periodic  motion  u  (t) 
with 


A  =  lim  A  ,  =  lim  A  . 

nl  n3 

n-*»  n-H» 


Au  =  lim  Au  ,  =  lim  Au  _ 

nl  .  n2 

n-w  n-*=° 


(1.21) 


In  order  to  demonstrate  this  intuitive  result ,  consider  the  situation  in 

Fig.  1.4.  The  loading  diagram  £(t)  and  the  value  of  k  are  the  same  as  in 

•  *  • 

Fig.  1.3,  but  initial  values  of  acceleration  uq  and  velocity  uq  have  been 
assumed  as  follows. 

The  values  of 


tn  ;  ti2  =  *11  +  2  :  t2i  =  tsii  +  T  ;  etc . 


(1.22) 


have  been  determined  first  from  the  condition 
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A* 

12 


A* 

13 


A* 

14 


A* 

n4 


(1.22  !  ) 


Then  u  and  u  have  been  assumed 
o  o 


so  that 


ii 

o 


u 

o 


-k 


(5  +  k)  dt 


'  o 


(1.23) 


(1.24) 


The  ensuing  motion  u(t)  ,  shown  in  Fig.  1.4,  is  periodic.  The  initial 

displacement  uq  is  arbitrary;  in  Fig.  L4c  it  has  been  chosen  such  that 

u  .  =  0  ,  so  that 

min 


0  $  u(t)  *  Au*  = 


+.A 

^n2 


nl 


u  dt 


(1.25) 


similar  to  (1.13). 

The  two  motions  in  Figs.  1.3  and  1.4  correspond  to  the  same  loading  and 

different  initial  conditions.  Therefore,  according  to  Martin’s  uniqueness 

theorem  [1966],  their  velocities  tend  to  each  other  as  time  increases.  In 

other  words,  the  actual  motion  u(t)  tends  to  a  periodic  motion  u  (t)  that 

s 

can  be  determined  as  in  Fig.  1.4  to  within  an  additive  constant,  u  .  ;  A 

mm  5 

g 

and  Au  in  Eqs.  (1.21)  are  given  by 


A*,  =  A* 
nl  n2 


(1.26) 


Au  =  Au* 


The  motion  u  (t)  will  be  indicated  as  the  steady  state  motion  (or  re¬ 
sponse);  the  actual  (limit)  value  of  u  .  as  the  final  residual  displacement 

mm  - 


1 R  * 
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g 

Note  that  the  rapidity  of  convergence  of  u(t)  to  u  (t)  ,  for  a  given 
load-time  diagram,  decreases  with  decreasing  k  .  In  the  limit,  if  k  s  0  , 
a  motion  with  a  periodic  (but  not  symmetric)  velocity  and  unbounded  displace¬ 
ments  is  activated  from  the  first  cycle*. 

After  a  sufficiently  large  number  of  cycles  the  large  strain  fatigue  can 
be  computed  on  the  periodic  motion  with  good  approximation 


u(n)  K  2nAu* 


(1.27) 


The  residual  displacement  after  n  cycles  u^n  is  given  by  (Fig.  1.4) 


n 


“to  =  Jx(4uil  -  Sui2>  *  n(4ull  -  412) 


(1.28) 


The  total  displacement  during  the  n**1  cycle  is  limited  by 


uR(n-l)  5  u(t)  5  uR(n-l)  +  %1 


(tnl  S  1  5  W 


(1.29) 


After  a  sufficiently  large  number  of  cycles,  Eq.  (1.29)  is  well  approximated  by 


u^  *  u  $  u^  +  Au*  (1.29') 

where  the  final  residual  displacement  uD  is  the  limit 

- —  ft 

<30 

UR  =  lim  uRn  =  L(Auil  "  Aui2}  (1-30) 

n-H»  1 

Finally,  note  that  if  in  Figs.  1.2,  1.3,  1.4  the  sign  of  the  5-diagrams  is 
changed,  the  signs  (and  the  signs  only)  of  the  u  and  u  diagrams  change  also. 

In  conclusion,  the  response  of  a  rigid-plastic  structure,  under  assumptions 
(1.05)  and  (1.06)jis  (or  tends  to)  a  periodic  motion  of  finite  amplitude  and 


One  cycle  of  such  a  response  is  illustrated  in  Fig.  4.7  at  the  end  of  the  re¬ 
port. 
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the  same  period  as  the  forcing  function  5(t)  .  At  each  instant  in  which  u 
becomes  zero,  there  is  a  discontinuity  6u  of  the  acceleration, 

-2k  s  6u  $  2k  (1.31) 

There  is  no  resonance  effect:  the  characteristics  of  the  response  are 
monotonic  functions  of  the  characteristics  of  the  load  input. 

It  is  possible  to  generalize  the  previous  results  to  periodic  but  not  sym¬ 
metric  loading  (only  the  first  of  Eqs.  (1.06)  holds)  and/or  non -symmetric  struc¬ 
tural  behavior  (k+  /  k  )  .  In  these  cases ,  a  permanent  displacement  is  added 
in  each  loading  cycle,  and  the  response  is  equal  (or  tends)  to  a  motion  sum  of 
a  periodic  plus  a  linear  term.  This  behavior  will  not  be  investigated  further, 
although  in  Section  4  a  sort  of  induced  asymmetry  of  structural  strength  will 


be  considered. 


-12- 


2.  "MULTI-MECHANISM"  RESPONSE  AND  "SHEAR  FRAMES”  SUBJECTED  TO  HORIZONTAL 
LOADS 


All  the  results  in  Section  1  were  obtained  in  the  assumption  that  the 
mechanism  (i.e.  the  "shape"  of  the  deformation  field)  did  not  change  in  time,  al¬ 
though  some  of  them  (and  in  particular  the  periodicity  of  the  steady  state  re¬ 
sponse  to  periodic  loading)  are  indeed  valid  also  if  this  assumption  is  lifted. 

In  the  latter  case,  graphical  operations  essentially  analogous  to  those 
illustrated  in  Figs.  1.2,  1.3  and  1.4  would  give  the  actual  and  the  periodic 
response.  However,  the  accelerations  corresponding  to  the  different  degrees 
of  freedom  are  given  by  relations  rather  more  complicated  than  Eqs.  (1.02)  and 
different  in  the  different  instances,  so  to  make  a  general  discussion  very  in¬ 
volved  if  not  impossible. 

In  fact,  the  departure  of  the  response  from  the  simple  single -me chan ism 
type  can  take  place  in  so  many  ways  that  a  unique  formulation  does  not  seem 
possible . 

To  illustrate  this  point,  it  is  sufficient  to  consider  a  uniform  simply- 
supported  beam  of  limit  moment  M^  ,  once  loaded  by  a  uniform  dynamic  load 
£(t)q  (Fig.  2.1a)  and  once  loaded  by  a  central  concentrated  force  £(t)P 
(Fig.  2. If).  In  the  first  case,  as  long  as  k  <  £  $  3k  ,  the  acceleration  and 
moment  diagrams  look  like  in  Figs.  2.1b  and  2.1c  respectively;  but  as  soon  as 
5  >  3k  ,  the  central  plastic  hinge  spreads  into  a  plastic  zone  of  finite  length 
Afc  ,  the  acceleration  (not  velocity!)  and  moment  diagrams  becoming  qualitatively 
as  shown  in  Figs.  2. Id  and  2.1e  respectively.  In  the  case  of  Fig.  2. If,  at 
first  the  acceleration  diagram  looks  again  like  in  Fig.  2.1g;  successive  stages 
(I,  II,  III)  of  the  bending  moment  diagram  are  illustrated  in  Fig.  2.1h.  At  a 
certain  value  of  £  ,  the  minimum  (max.  negative)  moment  M*  reaches  the  negative 
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limit  moment  -M  ,  and  new  plastic  hinges  are  formed:  further  acceleration 
and  moment  diagrams  look  like  in  Figs.  2.1i  and  2.1 j. 

Therefore,  instead  of  attempting  a  general  discussion  of  the  multi-degree- 
of -freedom  dynamic  rigid-plastic  behavior,  the  present  section  deals  with  a 
specific  structural  type  and  loading,  namely  rectangular  multi-storey  shear 
frames  subjected  to  horizontal  ground  motion  (as  shown  in  Fig.  2.2),  for  which 
the  symmetry  conditions  (1.05-06)  are  verified  and  particularly  relevant,  and 
moreover  the  motions  of  the  several  degrees  of  freedom  are  determined  by  re¬ 
markably  simple  relationships . 

The  following  discussion,  interesting  in  itself  as  the  solution  of  an  im¬ 
portant  structural  problem,  is  also  useful  as  an  illustration  of  the  techniques 
which  could  be  employed  for  other  analogous  problems. 

The  behavior  of  the  multi-storey  shear  frame  can  be  assimilated  to  that  of 
the  simple-supported  beam  in  Fig.  2. If,  in  that  a  second  degree  of  freedom  is 
activated  by  formation  of  new  plastic  zones  (or  plastic  hinges),  while  in  the 
beam  of  Fig.  2.1a  the  shape  of  the  deformation  fields  is  altered  by  the  spread¬ 
ing  of  the  plastic  zone.  Another  possibility,  substantially  different  from  both 
instances  in  Fig.  2.1,  will  be  illustrated  in  Section  3. 

2.1  "Shear  Frames";  Yielding  of  One  Storey 

We  define  by  shear  frame  a  rectangular  multi-storey  frame  such  that,  under 

horizontal  loads  acting  on  the  beams ,  the  plastic  hinges  form  in  the  top  and 

bottom  sections  of  the  columns.  The  rigid-plastic  deformations  of  such  a  frame 

are  fully  specified  by  N  independent  variables  u..  ,  where  u.  is  the  right- 

til  til 

ward  displacement  of  the  x  floor  with  respect  to  the  (i-l)x  and  N  the  num¬ 
ber  of  storeys  (Fig.  2.2). 


-14- 


The  shear  strength  of  the  i  storey 

J. 

I  .  M  . 

T  1  OD 

V  .  =  2  --—rr - 

oi  h. 

l 


(2.01) 


is  supposed  known  (and  equal  in  both  directions),  while  the  number  of  columns 
is  irrelevant,  and  can  indeed  vary  from  one  storey  to  another,  as  the  other 
quantities  on  the  right  side  of  Eq.  (2.01).  Only  plane  frames  are  studied  in 
this  section,  while  three-dimensional  shear  frames  will  form  the  object  of  Sec¬ 


tion  3. 


The  frame  will  be  assumed  to  be  subjected  to  a  horizontal  ground  movement 
whose  acceleration 


u  =  -£(t) 

g 


(2.02) 


is  a  known  function  of  time. 

The  mass  of  the  i**1  floor  is  itk  ,  and  the  mass  of  the  columns  is  negligi¬ 
ble  in  comparison.  Then  the  equivalent  loading  on  the  frame  is  constituted  by 
a  set  of  horizontal  forces  m.£  acting  at  each  floor  level,  and  the  total  in¬ 
ertia  forces  are* 


X 

7.  =  m.  (£  -  J.  vi.) 

i  i  p  D 


(2.03) 


The  total  shear  at  the  i  storey  is  given  by 


N  N  i 

vi  *  U  rz  -  L[V«  -  L  a  >: 


(2.04) 


Let  us  assume  that  at  a  certain  instant  the  1  storey,  but  none  above  it, 
is  deforming  plastically,  with  a  positive  velocity  u^  ,  i.e. 


Cf.  the  simplest  particular  case  in  Fig.  1.1. 
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y.  u.  =  y.  u.  for  Z  *  i 

*-*  "I  *1 


n  3 


v.  =  v  . 

i  oi 


(2.05) 


0  <  V  <  V  for  Z  >,  i 
oZ 


Because  of  (2.01)  and  (2.04),  Eqs.  (2.05)  implies 


V  . 

— °JL_  < 


o  Z 


N  N 

/  •  in  •  )  •  in  • 

n  3  p  3 


(Z  *  i) 


(2.06) 


This  relation,  which  is  independent  of  the  ground  acceleration  -£  and 
of  time  t  ,  shows  that  plastic  yielding  starts  in  the  p**1  storey,  defined  by 


i  °P 

k  =  =  mxn 

l.  m. 

P  H 


V  . 
oi 

N 

y.  m. 

i)  3 


(i  =  1,2,. ..N) 


(2.07) 


and  cannot  move  upwards  along  the  frame.  And  if  (2.07)  is  satisfied  by  p  =  1 
(i.e.  by  the  first  or  ground  storey),  as  it  is  very  likely  to  be  in  an  actual 
frame,  then  the  motion  of  the  (rigid-plastic)  frame  involves  only  one  degree  of 
freedom:  the  results  of  Section  1  apply  in  full  and  no  further  discussion  is 
necessary  for  such  a  frame . 

It  is  worthwhile  to  note  that  this  condition  can  be  verified  just  because, 
at  every  storey  except  the  first ,  there  are  infilling  walls  and  other  non-struc- 
tural  components  which  increase  the  actual  shear  resistance  . 

Although  the  considerations  above  have  been  formulated  with  reference  to  a 
specific  load  condition,  they  do  remain  valid  for  other  horizontal  loads  under 
not  too  restrictive  assumptions,  as  long  as  the  definition  of  "shear  frame"  holds. 
In  conclusion,  the  response  of  a  rigid-plastic  building  frame  to  horizontal  loads 
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is  very  likely  to  involve  yielding  of  the  first  storey  only:  use  of  this  con¬ 
clusion  as  a  greatly  simplifying  assumption  will  be  of  enormous  help  in  the 
study  of  other  aspects  of  rigid-plastic  frame  behavior  in  later  sections  of 
this  report . 


2 .2  "Shear  Frames"  with  More  Storeys  Yielding 

Consider  now  the  case  in  which  Eq.  (2.07)  is  satisfied  by  p  >  1  .  This 
"th 

p  storey  yields  as  soon  as  £(t)  crosses  either  of  the  lines  £  =  ±k^  ; 

again  it  can  be  assumed  without  lack  of  generality  that  £  =  +k^  occurs  first , 

th 

say  at  t  =  t  (Fig.  2.3).  For  t  1  t  ,  and  as  long  as  only  the  p  storey 

J.J.  21 

yields , 


u  =  C  -  jv-°P-  =  £  -  k 
p  N  p 

m. 

L1  j 

P 


^  =  0  (i  /  p) 


(2.08) 


These  equations  hold  as  long  as* 

Pr1 

V.  =  5  ^  +  <C 


N 

u  )  J .  m.  <  V  . 
P  pH  3 


oi 


(i  <  p) 


(2.09) 


Let 


V  -V 
oq  op 

— r - =  min 

pv 

)  .  m. 

qj  j 


V  .-V 
oi  op 


Py1 

).  m. 

P  3 


(i  =  1 . p-1) 


(2.10) 


and  note  that,  as  it  can  be  easily  proved,  Eq.  (2.07)  implies 


k  >  k 
9  P 


(2.11) 


Since  it  has  been  already  demonstrated  that  yielding  cannot  "move  upwards," 
it  is  only  necessary  to  consider  storeys  below  the 
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If  ?(t)  remains  always  smaller  than  k^  ,  no  other  storey  but  the  p 
is  ever  yielded,  and  the  response  of  the  frame  is  again  a  one-degree-of  freedom 
motion,  for  which  the  results  of  Section  1  apply. 

Suppose  now  that  a  certain  time,  say  t^^  (Fig.  2.3), 

5(t.-)  =  k  (2.12) 

15  q 

Then,  Eqs.  (2.08)  hold  only  in  the  interval 


*11  *  *  *  "15 


(2.12) 


while  for  t  >,  t^,.  it  can  be  easily  derived 


iii  =  0  (i  /  p,q) 


(2.13) 


i.e.  two  storeys  are  yielding:  in  other  words,  two  "simple  mechanisms"  are 
active  at  the  same  time ;  the  acceleration  of  the  first  mechanism  remains  con¬ 
stant,  while  that  of  the  second  is  graphically  immediate  (dot-shaded  in  Fig. 
2.3). 

It  would  be  quite  simple  to  derive  further  values  of  the  ground  accelera¬ 
tions  which  activate  further  mechanisms,  and  to  write  the  relevant  expressions 
for  the  storey  accelerations,  analogous  to  (2.13). 

Suppose  instead  that  {[(t)  satisfies  the  periodicity  conditions  (1.06)  and 
that  its  maximum  value  is  such  that  only  the  two  already  considered  mecha¬ 

nisms  can  be  activated. 
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The  motion  beyond  t  >  t^g  is  governed  by  Eqs.  (2.13)  as  long  as 


[*  .. 

u  =  u  dt  i  0 

q  q 

Jt« 


(2.14) 


i.e.  up  to  the  time  t^g  ,  defined  by 


(C  -  k  )  dt  *  0 

q 


(2.14') 


or  graphically,  in  the  example  of  Fig.  2.3, 


A2  =  A3  +  A4 


(2.15) 


At  t  -  t^g  ,  the  velocity  of  the  first  mechanism  is 


Up<t16)  * 


u  dt  =  A,  +  A„  : 
p  13’ 


(2.16) 


Eqs.  (2.08)  are  again  valid  beyond  t^g  ,  up  to  the  instant  t^2  defined  by 
Up(tl2 )  =  f  12  iip  dt  =  f  15  (C  -  kp)  dt  +  (  16  (kq  -  kp)  dt  +  f  12  (£-kp)  dt=  0 


or,  graphically,  by 


(2.17) 


A1  +  A3  A5 


(2.18) 


But  Eqs.  (2.15)  and  (2.18)  imply 


A1  +  A2  =  A4  +  A5 


(2.19) 


i.e.  the  time  t^2  ,  at  which  the  velocity  of  the  first  mechanism  up  becomes 
zero,  can  be  determined  ignoring  the  existence  of  the  second  mechanism.  (This 


result  holds  for  any  value  of  5(t12)  *  ^ 
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The  further  motion  beyond  t^  can  be  easily  constructed  by  similar 
arguments.  Analogous  to  Section  1,  if 


5(t12>  5  -*p 


(2.20) 


for  the  assumed  periodic -symmetric  5(t)  ,  the  negative  velocity  pulse  be¬ 
gins  at 


"13 


*  h  +  1 4  ti2 


(2.20') 


and  a  periodic  motion  is  thus  started,  with  the  same  period  as  the  forcing 
function  £(t)  ,  and  comprising  intervals  in  which  none,  one  and  two  mecha¬ 
nisms  are  active. 

If  Eq.  (2.20)  is  not  satisfied,  the  same  arguments  developed  in  Section 
1,  prove  that  the  response  tends  to  a  steady  state  periodic  motion  whose  ac¬ 
celeration  diagrams  can  be  most  easily  determined  taking  advantage  of  the 
property  (2.19),  as  shown  in  Fig .  2.4. 

Namely,  the  instants 


*11  •  *12  *  *?!  *  2  '  *11=  *!l  +  T  •' 


(2.21) 


at  which  u^  =  0  ,  are  first  determined  by  the  condition 


A*  +  A*  =  A*  +  A*  =  A*  +  A*  =  A*  +  A*  = 
1  2  4  5  6  7  9  10 


(2.22) 


The  second  mechanism,  in  each  cycle  of  the  steady  state  response,  is  activa¬ 
ted  at  t  =  t*  and  t  =  t*  ,  which  coincide  with  t*  and  t*  if 
n5  n7  *  nl  n2 


(2.23) 


as  in  the  example  in  Fig.  2.4.  This  second  mechanism  stops  at  the  instants 
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t*  and  t*  =  t*  +  £  (2.24) 

n6  n8  n6  2 

given  by 

A*  =  A*  +  A*  =  A*  r  A*  +  A*  (2.24') 

and  certainly  distinct  from  and  t(n+l)l  • 

The  derivation  of  the  response  of  each  relevant  mechanism,  discussed 
above,  leads  for  each  mechanism  to  the  same  conclusions  reached  in  Section  1. 

Finally,  note  also  that  the  symmetry  of  the  frame  behavior  with  respect 

to  the  sign  of  the  load  implies  only  that  positive  and  negative  limit  shear 

forces  (2.01)  are  equal.  This  condition  is  less  stringent  than  the  general 
one;  for  instance,  it  is  satisfied  if  the  moment-axial  load  yield  profiles 
of  the  columns  are  symmetric  with  respect  to  the  axial  load  axis  but  not  with 

respect  to  the  origin  (as  is  the  case  for  reinforced  concrete  columns);  and 

it  may  be  satisfied  even  if  the  individual  columns  have  different  positive 
and  negative  limit  moments,  provided  they  are  suitably  arranged. 

By  the  same  token,  we  can  add  constant  loads  on  the  structure  (e.g.  ver¬ 
tical  loads  on  the  columns),  and  thus  change  only  the  numerical  value  of 
and,  consequently,  of  V  .  (2.01),  provided  the  symmetry  of  the  latter  with 

respect  to  the  variable  (horizontal)  loads  is  not  affected.  (This  conclusion 
does  not  take  into  account  the  possible  instabilizing  effect  of  the  vertical 


loads,  which  will  be  considered  in  Section  4). 
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3.  THREE-DIMENSIONAL  "SHEAR  FRAMES”  SUBJECTED  TO  HORIZONTAL  LOADS 

As  it  has  been  pointed  out  in  the  previous  section,  in  most  actual  in¬ 
stances  of  rigid-plastic  shear  frames  only  the  first  storey  yields  under  hori¬ 
zontal  loads.  For  the  purpose  of  the  present  study,  the  frame  thus  reduces  to 
a  single  rigid  mass 

N 

m  =  Y.  m.  (3.01) 

l1  1 

supported  by  a  number  of  rigid-plastic  columns .  The  inertia  forces ,  being  ap¬ 
plied  at  a  finite  height  over  the  top  of  the  columns,  cause  differential  axial 
loads  in  the  columns:  if  the  effects  of  these  loads  are  neglected  (as  they 
have  been  implicitly  in  Section  2),  the  rigid  mass  can  be  further  reduced  to  a 
platform  of  vanishing  depth  (Fig.  3.1). 

This  Section  3  is  concerned  with  the  behavior  of  a  frame  thus  simplified, 
when  the  columns  are  not  all  complanar.  The  static  collapse  under  a  horizontal 
force  P  is  discussed  first  in  Subsection  3.1,  then  Subsection  3.2  deals  with 
the  dynamic  response  to  a  time-dependent  force  P  or  a  horizontal  ground  ac¬ 
celeration  u  .  Note  that  the  loads  P  and  u  have  been  underlined  (as  all 

-g  -  -g 

other  relevant  quantities  in  this  section  will  also  be),  in  order  to  stress 
their  two-component  vectorial  character. 

The  plastic  dynamic  behavior  of  three-dimensional  frames  of  a  similar  type 
has  never  been  examined  before,  to  the  writer's  knowledge,  although  it  seems  to 
be  a  realistic  model  of  buildings  in  which  architectural  or  functional  reasons 
prevent  bracing  in  any  vertical  plane,  at  least  at  the  first  (ground)  storey 
level.  Therefore,  the  problem  studied  herein,  besides  being  an  interesting  ex¬ 
ample  of  a  multi-mechanism  response  of  a  type  not  yet  considered,  is  also  of  di¬ 
rect  technical  relevance,  and  well  worth  some  further  attention.  Although  only 
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preliminary  results  are  presented,  there  seems  to  be  no  conceptual  difficulty 
in  their  application  and  generalization.  In  fact,  at  least  under  some  limita¬ 
tive  conditions,  the  rigid-plastic  hypothesis  already  allows  to  derive  some 
rather  general  qualitative  indications  and  to  define  some  properties  essential 
for  further  analyses  of  a  more  quantitative  nature. 


3 . 1  Static  Collapse 


The  static  collapse  of  the  frame  in  Fig.  3.1,  which  has  been  the  object  of 
recent  works  by  Harrison  [1968]  and  Wittrick  [1968],  can  be  described  as  an  ele¬ 
mentary  rotation  dp  of  the  rigid  platform  about  a  vertical  axis  c  .  In  the 
following,  it  will  be  assumed  (as  already  done  by  both  Harrison  and  Wittrick) 
that  the  torsional  strength  of  the  columns  is  negligible  in  comparison  with  their 
bending  strength  in  any  plane .  It  is  then  very  easy  to  find  the  collapse  value 

P  and  the  line  of  action,  f  ,  of  the  force  P  associated  to  each  axis  of  ro- 
-o 

tat ion  c  (although  sometimes  P^  is  not  defined  uniquely  for  a  given  c  ,  as 
it  will  be  seen). 

Let  C  and  (Fig.  3.2a)  be  the  traces  on  a  horizontal  plane  of  the 

axis  of  rotation  c  and  of  the  axis  of  the  column  respectively:  let  fur¬ 
ther  be  the  yield  profile  for  biaxial  bending  of  the  same  column*.  The 

axes  of  the  plastic  hinges  developing  in  the  top  and  bottom  cross-sections  of 
the  column  are  parallel  to  CB^  ;  the  corresponding  yield  moment  vectors  M  . 
are  obtained  from  the  normality  rule  (Fig.  3.2a);  the  shear  force  in  the  yield¬ 
ing  column  is  orthogonal  to  the  moment  vectors  and  of  magnitude 


M  . 

V  .  =  2  £1 

o]  h 


(3.02) 


for  a  column  of  length  h  built-in  at  both  ends. 


Yield  profiles  of  structurally  significant  sections,  including  the  effect  of 
axial  load,  have  been  recently  presented  by  S.  Santathadapom  and  W.  F.  Chen 
[1968]. 
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Alternatively  (Fig.  3.2b),  the  moment  yield  profile  can  be  rotated  by  90° 
and  divided  by  h/2  to  obtain  the  shear  yield  profile  of  the  j  column.* 

The  normality  rule  can  then  be  applied  directly  to  obtain  the  shear  yield  force 


from  the  elementary  displacement 


du.  =  dp  x  CB. 
-D  -  —3 


(3.03) 


This  alternate  procedure  is  slightly  more  general,  because  it  includes  the  cases 
in  which  the  shear  yield  profile  is  obtained  less  immediately  from  the  moment 
profile  (e.g.,  when  the  bottom  section  of  the  column  is  pinned  in  one  direction 
only,  or  when  the  top  and  bottom  sections  of  the  column  are  different  from  each 
other) . 

In  the  more  general  situation,  the  vectorial  sum  of  the  Vq.  gives  uniquely 
the  collapse  force  associated  with  C 


P 

-o 


=  l.  V  . 
L1  "03 


(3.04) 


Indeterminacies  arise  whenever  a  shear  yield  profile  has  a  straight  segment  nor¬ 
mal  to  diu  ,  and  in  the  particular  case  of  C  coinciding  with  one  of  the  , 

say  with  B 
J  n 

th 

In  the  last  case,  the  n  term  in  the  sum  (3.04)  is  only  required  not  to 
violate  and  is  otherwise  arbitrary.  Wittrick  [1968]  has  shown  that  this 

situation  arises  when  the  foot  of  the  perpendicular  from  B^  to  f  falls  within 
a  plane  domain,  called  safe  zone  of  the  n ^  column  and  given  by  a  very  simple  geo¬ 
metrical  transformation  of 

Vn 


The  shear  yield  profile  thus  defined  is  not  connected  with  a  ,rshear  fail¬ 
ure”  of  the  column,  but  with  the  ”shear  force”  which,  at  the  collapse  limit, 
equilibrates  the  moments  in  the  plastic  hinges  formed  at  the  ends  of  the  col¬ 
umn. 


On  the  other  hand,  if  C  is  at  infinity,  i.e.  if  the  collapse  motion  of 
the  platform  is  a  translation,  all  du^  are  equal,  say  to  du  :  however,  the 
associated  are  in  general  different  and  not  parallel  to  each  other.  By 

varying  the  direction  of  the  translation  du  ,  the  value  and  the  line  of  action 
of  Pq  (3.04)  are  changed.  If  the  Pq  associated  with  all  possible  transla¬ 
tions  are  drawn  with  origin  in  the  same  point,  their  end  points  describe  a 
closed  convex  curve ,  symmetric  with  respect  to  the  origin ,  that  can  be  called 
total  shear  yield  profile  (or  more  simply  yield  profile)  of  the  frame  $p  :  the 
normality  rule  between  4>p  and  du  holds. 

On  the  contrary,  it  does  not  appear  possible  to  give  an  evident  geometrical 
condition  or  property  satisfied  by  the  lines  of  action  f  of  Pq  corresponding 
to  translational  collapse.  This  is  shown  by  the  comparatively  very  simple  ex¬ 
ample  in  Fig.  3.3,  pertaining  to  a  platform  supported  on  three  I-section  columns 
,  b2  ’  ®3  w^ose  yield  profiles  are  ,  $V2  ,  $V3  respectively:  each 

translation  vector  du^  and  the  corresponding  line  f^  are  labeled  by  the  same 
suffix . 


Fortunately,  in  a  number  of  instances  of  practical  interest  a  point  ex¬ 
ists  such  that,  when  f  passes  through  with  any  inclination,  the  collapse 

motion  is  a  translation.  Such  a  point  can  be  well  defined  as  the  plastic 

shear  center  of  the  frame.  (Note  however  that  is  defined  by  a  purely 

sufficient  condition,  which  does  not  exclude  that  a  Pq  with  line  f  not  throigh 
(1!^  may  also  cause  translational  collapse.)  Examples  of  are  as  follows: 


a)  if  the  plan  of  the  columns  has  two  axes  of  symmetry,  co¬ 
incides  with  their  intersection; 
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oriented,  any  du  is  (or  can  be)  associated  with  yield  forces 

V  .  parallel  to  each  other  and  in  ratios  s.  :  C  is  there- 
-03  3  P 

fore  the  centroid  of  masses  proportional  to  s^  concentrated 


in  B. 


circles . ) 


(This  example  is  always  verified  when  the  ty.  are 


c)  if  the  are  rectangles  defined  by 


V  .  *  V  . 

X3  0x3 


(3.05) 


V  .  <  V  . 

Y3  oy3 


C  is  given  by  the  intersection  of  the  resultants  of  V  .  and 
P  0x3 

V  . 

oy3  • 

3.2)  Dynamic  Behavior 

The  behavior  of  the  frame  in  Fig.  3.1,  subjected  to  a  time -dependent  force 
applied  in  the  centroid  CL 


P  =  m  C(t) 


(3.06a) 


or  to  a  horizontal  ground  acceleration 


u  =  -£(t)  (3.06b) 

-g 

will  now  be  examined.  For  the  sake  of  simplicity,  in  the  following  discussion 
it  will  be  assumed  (i)  that  the  magnitude  and  the  sign,  but  not  the  direction 
of  £  ,  vary  in  time,  and  (ii)  that  the  plastic  shear  center  exists  and 

coincides  with  the  centroid  of  the  platform  C„  (as  it  does  in  example  a)  at 

b 

the  end  of  Section  3.1).  The  frame  is  thus  loaded  at  any  instant  by  a  total  re¬ 
sultant  force  applied  in  the  plastic  shear  center 


F  =  m(5  -  u) 


(3.07) 
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where  u  is  the  acceleration  of  the  platform  relative  to  the  base  plane  (Fig. 
3.4a).  The  motion  is  a  translation  whose  instantaneous  velocity  u  satisfies 
the  normality  rule  with  respect  to  the  yield  profile  of  the  frame  $p  .  In 
analogy  with  Sections  1  and  2,  it  is  convenient  to  divide  the  yield  profile  by 
m  ;  the  profile  thus  obtained  will  be  denoted  by  $  and  its  equation  symboli¬ 
cally  indicated  by 

<f>(k  »  k  )  =  0  (3.08) 

x  y 

(Fig.  3.4b).  Then,  letting 


ox  =  |  =  5(?x,  y 

••  M.  •*  ••  - 

LX  =  u  =  u(ux,  uy)  , 


(3.09) 


the  conditions  governing  the  motion,  analogous  to  Eqs.  (1.02),  can  be  stated  as 
follows : 


•  •• 

a)  when  u  i  0  and/or  u  f  0  ,  the  total  loading  point  L  is 
on  the  yield  curve  $  ,  and  u  is  directed  as  the  outward 
normal  to  t  in  L  (Fig.  3.4b); 


*(kx>  ky>  =  -  V  V  -  y  *  0 


(!!_> 

v3k  ;| 


u  + 
xo 


(£  -k  )  dt 
x 


L  Ux 


‘SM 


Ci  rt 

y  •  1 

y  u  + 

y° 


(^  -k  )  dt 

y  y 


(3.10) 


where  t  is  the  initial  instant  of  motion  and  u  ,  u  the  components  of 
o  xo  *  yo  r  . 

a  possible  initial  velocity  uq  .  If  L  coincides  with  a  pointed  vertex  of  $ 
the  usual  well  known  generalizations  of  Eqs.  (3.10)  apply. 


9 
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b) 


when  u 


0  and  X  is  not  external  to  <J>  , 


♦(5X.  5y)  ^  0 

u  =  0  (X  =  L) 


(3.11) 


All  the  information  on  the  motion  might  be  obtained  from  two  plots ,  analo¬ 
gous  to  those  used  in  Sections  1  and  2,  namely 


a) 

U  =  £ 

-  k 

vs .  t 

X  X 

X 

b) 

u  =  £ 

-  k 

vs.  t 

y  y 

y 

where  £x(t)  and  £y(t)  are  given  and  proportional  to  each  other  by  assump¬ 
tion,  k^(t)  and  ky ( t )  are  such  that  either  (3.10)  or  (3.11)  are  verified. 
However,  the  variation  of  k^  and  ky  with  t  is  not  known  a  priori  and  the 
diagrams  (3.12)  can  be  in  general  obtained  only  by  some  form  of  successive  ap¬ 
proximation  procedure. 

For  the  following  discussion,  it  is  useful  to  define  a  principal  direction 
of  the  yield  profile  as  the  direction  of  a  straight  line  ru  through  the  ori¬ 
gin,  orthogonal  to  the  yield  profile  in  the  points  of  intersection  I,.  and  I| 
(e.g.  ,  n2  in  Fig.  3.5a).  If  1^  and  I!  are  vertices  of  the  yield  pro¬ 

file,  the  definition  is  generalized  by  requiring  n  to  be  comprised  in  the  fan 
of  normals  to  the  yield  profile  in  I.  or  I!  (e.g.  n.  ,  ri  ,  P-  in  Fig. 

11  £.  O  O 

3.5b).  Consider  now  a  straight  line  ru  slightly  rotated  around  0  with  re¬ 
spect  to  n.  :  the  tangential  component  8^  of  an  outward  vector  directed  along 
can  either  point  towards  the  original  position  tk  ,  or  away  from  it  (e.g., 
respectively  0^  and  8^  in  Fig.  3.5a).  If  the  first  condition  is  verified  on 
both  sides  of  tk  ,  it  will  be  said  that  the  principal  direction  tk  is  stable; 
viceversa,  n.  is  said  to  be  unstable  if  the  second  condition  is  verified,  at 
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least  on  one  side .  The  physical  reason  for  the  use  of  these  terms  will  soon 
become  apparent:  examples  of  unstable  principal  directions  are  n ^  in  Fig. 
3.5a,  TI2  »  n3  »  n5  in  fig.  3.5b.  Because  of  the  symmetry  of  ♦  with  re¬ 
spect  to  the  origin,  at  least  two  principal  directions  always  exist,  of  which 
one  is  stable  and  the  other  unstable. 


If  the  loading  vector  £  is  directed  along  a  principal  direction  n  , 

•  «  • 

also  the  acceleration  u  and  the  velocity  u  are  at  all  times  parallel  to  n  ; 
the  total  loading  point  L  lies  in  the  segment  of  n  not  external  to  the 
yield  profile.  In  other  words,  the  response  involves  only  one  degree  of  free¬ 
dom  ,  the  two  diagrams  (3.12)  are  proportional  to  each  other,  and  the  whole  dis¬ 
cussion  in  Section  1  is  fully  applicable.* 

On  the  contrary,  if  |  is  not  parallel  to  a  principal  direction,  the  ac- 

••  •  # 
celeration  u  and  the  velocity  u  change  direction  with  time,  and  the  motion 

is  a  two -degree -of -freedom  translation.  If  the  vectors  u  at  successive  in¬ 


stants  are  reported  from  a  common  origin  0^  , 


0  U  =  u  =  u  +  udt 
v  -  -o 


(3.13) 


(where  a  vector  integral  appears),  their  end  point  U  describes  a  trajectory  v  , 
whose  tangent  in  any  instant  is  parallel  to  the  corresponding  u  . 

As  a  first  example,  consider  the  response  to  an  impulsd+corresponding  to  an 

1  1 

initial  velocity  u  applied  at  time  t  =  0  (Fig.  3.6a).**  For  any  t  >  0  , 


If  §  is  parallel  to  n  at  every  storey,  the  developments  in  Section  2  are 
valid  for  a  multi-storey  three  dimensional  shear  frame  as  well.  Otherwise, 
the  latter  structure  must  be  analyzed  by  generalizations  of  the  procedures 
outlined  in  this  section. 

An  impulsive  loading  problem  of  a  similar  type ,  applied  to  a  lumped-mass 
cantilever  beam,  has  been  the  subject  of  a  very  recent  paper  by  Frick  and 
Martin  [1968]. 
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5=0,  i.e.  X  =  0  (3.14a) 

and ,  in  the  interval  0  <  t  <  t^. 

u  =  LO  (3.14b) 

The  initial  position  of  L  is  determined  by  the  requirement  of  normality 

between  $  and  u  and  gives  the  initial  tangent  to  the  v-trajectory  (Fig. 

-o 

3.5b).  From  the  qualitative  construction  of  v  reported  in  Fig.  3.6b,  it  is 
evident  that  the  final  tangent  to  v  is  parallel  to  a  stable  principal  direc¬ 
tion.  In  other  words,  after  the  impulse,  the  total  loading  point  L  moves 
along  the  yield  curve  in  the  direction  determined  by  the  relevant  tangential 
component  of  uq  ;  when  the  motion  stops ,  L  has  reached  but  not  passed  the 
first  stable  point  of  intersection  1^  . 

It  is  thus  evident  that,  if  the  direction  of  the  impulse  is  slightly  ro¬ 
tated  with  respect  to  a  stable  principal  direction,  the  response  varies  only 
slightly;  while  even  a  very  little  deviation  of  the  load  vector  from  an  unsta¬ 
ble  principal  direction  causes  a  large  difference  at  least  in  some  aspects  of 
the  response  (like  the  movement  of  L  along  $).  In  other  words,  the  one- 
degree  -of  -freedom  response  is  fully  reliable  only  when  the  loading  vector  is 
directed  along  a  stable  principal  direction.  This  property,  which  appears  to 
hold  for  any  loading  history,  justifies  the  terms  stable  and  unstable ;  however, 
it  is  not  yet  clear  in  general  terms  how  the  variation  in  the  movement  of  L 
affects  the  structural  response.* 

The  motion  of  L  along  $>  need  not  be  regular  or  continuous.  Take  for 
instance  the  most  simple  case  of  impulsive  loading  illustrated  in  Fig.  3.7. 
Initially,  L  remains  in  Lo  for  a  finite  time  interval,  while 

ft 

Note  also  that  most  of  the  considerations  just  developed  with  regard  to  the 
load  5-u ,  yield  profile  and  velocity  u  of  this  particular  structural  model, 
apply  'also  to  the  stress  vector,  yield  surface  and  strain  vector  in  each  point 
of  a  rigid -perfectly  plastic  three-dimensional  continuum. 
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u(ii  ,  u  )  =  L  0  (-k,  ,  -k„)  =  const 
-  y*  x  o  1’  2 


(0  <  t  <  tx) 


(3.15) 


At  t  =  t^  ,  the  velocity  vector  u  becomes  parallel  to  x  ;  L  jumps  to  , 

•  •• 

then  stays  there  while  both  u  and  u  remain  horizontal, 


u(u  ,  u  )  =  L,0  (0,  -k_)  =  const 

_  y  7  J_  7  2 


(3.15* ) 


(t1  <  t  <  tf) 


The  motion  finally  stops  at  t  =  t^.  .  The  v-trajectory  is  constituted  by  two 
straight  segments  (Fig.  3.7b);  the  velocity-time  diagrams  are  immediate  (Fig. 
3.7c,  d) . 

As  a  further  example,  consider  a  frame  whose  strength  in  one  particular 
direction,  say  x  ,  is  much  larger  than  both  the  strength  in  the  orthogonal  di¬ 
rection  y  ,  and  the  maximum  x-component  of  external  load,  ^xm  •  't^^s 

case,  the  yield  profile  can  be  approximated  by  two  parallel  straight  lines  (Fig. 
3.8a) 


k  =  ±k. 

y  1 


(3.16) 


and,  although  both  X  and  L  vary  with  time,  u  and  u  are  parallel  to  y 
for  any  loading  history,  i.e. 


u=£-k=0;  u=0 

XXX  X 


(3.17) 


The  motion  can  be  studied  as  in  Section  1  for  the  non-vanishing  component  (Fig. 
3.8b) 


U  =  5  -k  =  e  +  k, 

y  y  y  y  i 


(3.17') 
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Finally,  if  the  yield  profile  of  the  frame  is  a  rectangle  with  sides  par¬ 


allel  to  x  and  y  ,  as  in  example  a)  at  end  of  Section  3.2,  the  two  compon¬ 


ent  diagrams  (3.12)  can  be  studied  independently  of  each  other  following  the 


procedures  outlined  in  Section  1.  In  particular, 


C +k.  when  u  >  0 

I  1  y 

k  =  <  .  (3.18a) 

y  1  -k,  when  u  <  0 

'-l  y 


C  +k_  when  u  >  0 

12  x 

kx  = <  #  (3.18b) 

1  -k.  when  u  <  0 

v.  2  x 

Both  x  and  y  are  stable  principal  direction  of  the  yield  profile,  its  diag¬ 
onals  are  unstable  principal  directions . 

If  the  maximum  x-component  of  the  load,  ,  (Fig.  3.9)  is  not  larger 

than  k^  ,  the  response  is  the  same  one-degree-of-freedom  just  described  with 
reference  to  Fig.  3.8. 

On  the  contrary,  Figs.  3.9  and  3.10  describe  qualitatively  the  steady- 
state  velocity  under  a  periodic -symmetric  loading  history,  such  that  both  com¬ 
ponents  fall  in  Case  II  defined  in  Section  1.2.  (Note  that,  by  an  arbitrary 
choice  of  scale,  the  ([^(t)  and  £y(t)  diagrams  have  been  drawn  equal.)  In¬ 
spection  of  the  velocity  diagrams  shows  that  u  never  vanishes:  the  normality 
rule  requires  the  total  loading  point  L  to  jump  from  one  vertex  to  the  other 
around  the  yield  profile  (from  to  L to  Lg  to  to  ,  etc.),  re¬ 

maining  in  each  vertex  for  a  finite  time  interval,  indicated  at  the  bottom  of 
Fig.  3.10.  The  curve  v  ,  already  defined,  is  shown  in  Fig.  3.9b:  it  is  (in 
the  steady-state  motion)  a  closed  curve,  symmetric  with  respect  to  the  origin 

0  . 
v 
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The  motion  of  L  around  the  yield  profile  is  slightly  more  complicated 
if  one  of  the  component  motions  (or  both)  falls  in  Case  I  (Section  1.1),  as  in 
the  example  qualitatively  illustrated  in  Figs.  3.11  and  3.12.  Namely,  L 
jumps  from  not  directly  to  ,  but  to  an  intermediate  point  L ^  »  then 

it  moves  gradually  towards  L_  .  in  the  interval  (t  ,  t  .),  while  both  u  and 

c.  Xz  Xo  X 

iix  remain  zero.  The  same  happens  between  and  . 

In  all  described  examples  (Figs.  3.7,  3.9,  3.11),  the  jumps  of  the  total 
loading  point  L  along  the  contour  of  the  yield  profile  $  (i.e.  discontinui¬ 

ties  in  the  acceleration  vector  ii  =  LX)  are  connected  with  the  existence  of 
straight-line  segments  in  $  ,  and  the  ensuing  discontinuities  in  the  points 
associated  with  a  continuously  varying  velocity  vector. 

Results  essentially  similar  to  those  derived  in  the  above  discussion  can 
be  also  obtained  with  more  general,  curvilinear  yield  profiles,  although  the 
actual  computation  of  the  whole  motion  would  probably  require  a  rather  compli¬ 
cated  successive  approximation  procedure.  It  seems  however  possible  to  state 
that  the  steady  state  response  of  the  three-dimensional  frame  under  considera¬ 
tion,  when  subjected  to  a  horizontal  periodic-symmetric  load,  is  a  periodic 
translation  with  the  same  period  as  the  loading,  of  which  three  subcases  can 
be  distinguished: 

a)  if  the  load  is  directed  along  a  principal  direction  of  the 
yield  profile,  the  response  is  a  one-degree -of- freedom  trans¬ 
lation,  parallel  to  the  same  direction; 

b)  if  the  direction  of  the  load  is  sufficiently  close  to  a  stable 
principal  direction,  the  total  loading  point  L  moves  in  two 
distinct  intervals  of  the  contour  of  ♦  ,  and  the  velocity 
vanishes  twice  during  each  cycle ;  the  v-curve  has  an  8-shape ; 
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c)  if  neither  a)  nor  b)  are  verified,  the  velocity  never  van¬ 
ishes  ;  L  moves  around  the  yield  profile ,  possibly  with 
some  jumps .  The  signs  of  the  rotations  of  the  total  load¬ 
ing  point  L  around  the  yield  profile  ♦  and  of  the  ve¬ 
locity  point  U  around  the  v-trajectory  appear  to  be  both 
determined  by  the  sign  of  the  angle  between  the  direction 
of  £  and  the  outward  normals  to  $  in  the  points  of  in¬ 
tersections:  e.g.,  it  is  clockwise  when  the  loading  vec¬ 
tor  is  directed  as  £^  in  Fig.  3.13,  anticlockwise  when 
the  load  is  directed  as  ^  • 

Note  that  there  is  no  continuity  between  the  motions  of  L  associated 
with  loads  directed  along  an  unstable  principal  direction  and  along  a  slightly 
rotated  direction. 

In  conclusion,  the  considerations  developed  in  this  section,  although 
still  very  far  from  practical  applicability,  have  yielded  some  general  indi¬ 
cations  on  a  two-degree-of -freedom  response  to  dynamic  load,  in  which  the 
change  of  mechanism  is  due  to  the  interaction  between  two  generalized  stress 
components  and  not  to  the  formation  or  spreading  of  plastic  zones.  The  dis¬ 
tinction  between  stable  and  unstable  principal  directions  can  be  of  great  in¬ 
terest  in  actual  problems . 

It  appears  that  successive  approximation  procedure  could  be  set  up  with¬ 
out  excessive  difficulties,  for  the  computation  of  the  motion  of  the  specific 
type  of  frames  so  far  examined,  in  which  it  will  be  recalled  that  by 

assumption . 

If  0^  does  not  exist  or  does  not  coincide  with  ,  the  response  is 
much  more  complicated,  in  that  it  involves  a  third  degree  of  freedom,  namely 
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a  rotation  about  a  vertical  axis.  Although  this  problem  has  not  been  investi¬ 
gated  at  all,  no  particular  obstacle  seems  likely  to  arise  in  the  extension  of 
the  arguments  developed  above,  with  the  obvious  exception  of  greatly  increased 
computational  difficulties. 
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4.  EFFECTS  OF  GRAVITY  ON  THE  DYNAMIC  BEHAVIOR  OF  SHEAR  FRAMES 

In  a  frame  such  as  those  examined  in  Section  2,  at  each  storey  level  the 
columns  are  subjected  to  compressive  axial  loads  due  to  the  weight  of  the  upper 
floors,  and  the  limit  moments  are  consequently  modified.  As  already 

pointed  out  at  the  end  of  Section  2,  this  effect  can  be  easily  accounted  for 
in  the  dynamic  analysis  developed  there,  just  by  substituting  the  modified 
limit  moments  in  the  expressions  of  the  frame  shear  strengths  (2.01)  at  each 
storey  level  (provided,  of  course,  that  the  shear  strengths  are  not  further 
significantly  modified  by  the  additional  column  axial  loads  due  to  the  hori¬ 
zontal  loads ) . 

The  present  section  will  rather  be  concerned  with  the  combined  instabiliz- 
ing  effect  of  gravity  loads  W  and  horizontal  displacements  u  in  a  plane 
shear  frame,  whose  dynamic  response  involves  yielding  of  the  lowest  storey  only 
and  which  therefore  can  be  diagrammatically  represented  as  in  Fig.  4.1.  (The 
number  of  columns,  as  pointed  out  already  in  Section  2,  is  irrelevant).  The 
following  discussion  could  easily  be  extended  to  cover  analogous  instability 
effects  in  the  single-mechanism  dynamic  responses  of  other  structures. 

With  reference  to  Fig.  4.1  as  compared  to  Fig.  1.1,  Eos.  (1.02)  are  modi¬ 
fied  into 

u  =  K  -  k  +  Gu  when  u  >  0 

u  =  K  +  k  +  Gu  when  u  <  0  (4.01) 

u  =  0  when  u  =  0  and  -(k  +  Gu)  £  £  S  (k  -  Gu) 

where 

e  =  -  or  e  =  -ii  (4.02) 

m  g 

measures  the  horizontally  applied  load,  and 
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(4.03) 


EW.  _  g 
mh  ~  h 


in  turn,  g  is  the  acceleration  of  gravity.*  Other  expressions  for  k  and 
G  would  be  appropriate  for  other  structural  types . 

In  the  ?  vs.  t  diagrams,  already  several  times  used  in  the  present 
report ,  the  effect  of  gravity  loads  is  a  vertical  translation  Gu  of  the 
strength  interval  (+k,  -k),  downwards  when  u  is  positive,  upwards  when  u 
is  negative  (Fig.  4.2). 


4.1  Impulsive  Load 

Consider  first  the  response  to  an  impulsive  load,  i.e.  to  an  initial  posi- 
tive  velocity  uq  .  Eqs.  (4.01)  give 


u  -  Gu  =  -k 


(4.03) 


for  any  t  >  0  ,  as  long  as  u  £  0  .  The  solution  of  Eq.  (4.03),  with  the  ini¬ 
tial  conditions 


u(0)  =  0  ;  u(0)  =  uq  >  0 


(4.04) 


yields 


V  u  _ 

7T  (1  -  cosh  *^Gt )  +  —  sinh  /Gt 
G  /G 

u  cosh  /§t  -  —  sinh  /Gt 
°  /G 

-k  cosh  /Gt  +  u  i^G  sinh  v^Gt 
o 


(4.05) 


In  a  building  frame,  h  is  of  order  10-15  ft;  hence 


G  -  2-3  sec  approx. 
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which  easily  reduce  to  (1.03)  when  G  0  . 

The  response  velocity  u  vanishes  at  a  time  t  =  t^  ,  given  by 

u  /G 

tanh  /Gt^  =  — (4.06) 

Since 

0  $  tanh  x  <  1  when  0  $  x  <  “ 


Eq.  (4.06)  yields  a  real  positive  total  response  time  t^  only  if 


u  /G 

<  1  (4.08) 


If  (4.08)  is  not  verified,  the  velocity  successive  to  the  impulse  never  van¬ 
ishes,  i.e.  the  frame  never  comes  to  a  stop  and  deforms  up  to  complete  col¬ 
lapse  .  * 

This  result  can  be  presented  in  several  alternative  ways .  For  example , 
introducing  (4.03),  disequality  (4.08)  becomes 


(m  u  )2  h ^ — -  <  1  (4.09) 

°  (EM)2 

o 

and  yields a  limit  relationship  between  the  sum  of  the  column  flexural  strengths 
(possibly,  modified  by  the  axial  loads),  the  storey  height  h  ,  and  the  impulse 
momentum  m  uq  that  a  rigid-plastic  shear  frame  can  absorb  without  complete 
collapse . 

u^ 

The  variation  of  velocity  with  time  for  several  values  of  — ,  and  the 
u  /g 

relationship  between  — and  the  total  response  time,  given  respectively  by 


* 

Remember,  however,  that  the  previous  equations  remain  quantitatively  valid 
only  as  long  as  the  displacement  u  is  small  in  comparison  with  h  . 
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the  second  Eq.  (4.05)  and  by  Eq.  (4.06),  are  plotted  with  solid  lines  in  Fig. 

4.3;  the  dotted  lines  represent  the  relationships  valid  when  the  instabilizing 
effect  is  neglected.  The  difference  of  behavior  when  (4.08)  holds  and  when  it 
does  not  is  very  evident:  in  the  latter  case,  the  velocity  decreases  at  first, 
reaches  a  minimum  and  then  increases  monotonically  up  to  the  ruin  of  the  frame.* 

4.2  Periodic  Load 

The  rest  of  the  present  section  deals  with  the  response  given  by  Eqs.  (4.01) 
when  the  load  function  £  *  £(t)  is  periodic-symmetric,  as  already  defined  by 
Eqs.  (1.06).  Comparing  Fig.  4.2  with  the  consideration  developed  in  Section  1, 
it  is  immediate  that  the  response  u(t)  ,  although  involving  successive  plastic 
deformations  of  alternate  signs,  cannot  be  periodic:  in  fact,  as  soon  as  u  t  0  , 
an  asymmetry  of  the  effective  strength  tends  to  increase  the  residual  displace¬ 
ments  of  a  defined  sign ,  which  build  up  during  the  successive  loading  cycles ,  un¬ 
til  the  critical  displacement 

u  =  u(t  )  =  ±  ~  (4.10) 

C  C  b 

is  reached ,  at  which  point  the  strength  of  the  frame ,  with  respect  to  horizontal 
loads  of  either  positive  or  negative  sign,  becomes  zero.** 

The  critical  time  (or  time  to  collapse)  tfi  can  by  physically  defined  by 
the  ability  acquired  by  the  gravity  loads,  to  cause  by  themselves  collapse  of 


Similar  results  were  obtained  by  Lee  and  Martin  for  a  rigid-plastic  column 
subjected  to  a  particular  type  of  transverse  impulse  [1966],  but  at  the  ex¬ 
pense  of  rather  cumbersome  algebra  that  made  the  ultimate  results  much  less 
immediately  evident  than  the  present  ones. 

Strain-hardening  with  Bauschinger  effect  would  also  cause  a  translation  of 
the  strength  interval  (-k,  +k),  but  in  a  stabilizing  direction;  the  response 
would  tend  to  a  periodic  response,  in  which  even  the  displacements  are  sym¬ 
metric.  Thus,  if  both  strain -hardening  and  gravity  loads  are  significant, 
their  effects  are  opposite  to  each  other,  and  either  can  prevail.  [August!, 
1969]. 
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the  frame.  In  other  words,  if  a  periodic  loading  £(t)  is  applied  for  a  fi¬ 
nite  number  of  cycles  n  ,  the  frame  collapses  (because  of  the  constant  axial 
loads)  if,  and  only  if, 

nT  >  t  (4.11) 

c 

or  otherwise ,  if 

n  >  n  (4.11* ) 

c 

where  the  critical  number  of  cycles  nc  is  defined  by 

n  T  $  t  <  (n  +  1)T  (4.12) 

c  c  c 

There  is  no  conceptual  difficulty  in  the  numerical  or  analytical  step-by- 
step  integration  of  Eqs.  (4.01)  for  any  given  loading  program,  and  in  the  di¬ 
rect  determination  of  the  critical  time  or  critical  cycle  nc  ,  which  will 

of  course  depend  on  the  characteristics  of  the  frame  and  on  the  shape  of  the 
loading  function.  Essentially  this  approach  has  been  used  in  recent  papers  by 
Jenning  and  Husid  [1968]  (where  an  elasto-plastic  frame  under  an  earthquake 
load  is  examined,  following  earlier  works  by  Housner  [1959,  1960]  who  appears 
to  have  been  the  first  to  discuss  the  phenomenon  of  collapse  of  plastic  struc¬ 
tural  frames  because  of  gravity  loads)  and  by  Ballio  [1968]  (who  found  an  anal¬ 
ogous  behavior  for  a  one -degree -of -freedom  elasto-plastic  column  model  and  a 
simple  periodic  loading).  The  numerical  process  would  be  greatly  simplified  in 
any  case  by  the  rigid-plastic  assumption:  the  response  could  be  determined, 
cycle  by  cycle ,  by  a  successive  approximation  procedure ,  as  outlined  below  with 
reference  to  Fig.  4.4. 

Suppose  that  the  displacement  un_1  [hence,  the  effective  strengths  (k  - 
-Gun_^)  and  (-k  -  Gun_^)]  at  the  beginning  of  the  n**1  cycle  are  known.  Assume, 
as  a  first  approximation,  a  straight  line  r  for  the  further  variation  of  Gu  . 
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and  perform  the  integrations  described  in  Section  1  from  t  up  to  a  tenta¬ 
tive  t^  :  introduce  the  diagram  of  u  ,  thus  obtained,  in  the  acceleration 
diagram  (Fig.  4.4a),  and  repeat  the  operations  until  the  resulting  t  ^  and 

Au  ,  are  stabilized.  The  same  procedure  is  then  used  to  obtain  t  and  Au_, 
nl  n4  n2 

It  can  be  noted  that ,  if 

5  <  2k  (4.13) 

m 

the  negative-velocity  phases  disapper  (i.e.,  Au  j  =  0  )  before  the  critical 
displacement  uc  is  reached. 

4.3  Direct  Determination  of  the  Critical  Cycle  nc  Under  Periodic  Loading 


For  any  practical  structure ,  nc  can  be  expected  to  be  rather  large :  in 
fact,  if  a  structure  is  at  all  intended  to  undergo  periodic  loading,  it  is  cer¬ 
tainly  required  to  be  able  to  withstand  more  than  a  few  cycles.  In  this  case 
the  rigid-plastic  assumption  appears  to  become  really  convenient,  in  that  it 
will  allow  the  direct  determination  of  bounds  to  nc  without  following  the 
whole  response  history,  as  outlined  below.  It  is  still  an  open  question  whether 
the  critical  times  thus  calculated  are  realistic,  although  rough,  approximations 
of  the  critical  time  of  the  actual,  elastic-plastic  frame. 

The  case  corresponding  to  Gase  I  of  Section  1.1  is  discussed  first. 

It  is  self-evident  from  Fig.  4.4  that  the  net  variation  of  u  in  each 

cycle 


Vl  =  4unl  -  lun2 


(4.14) 


increases  with  n  :  hence 


(4.15) 
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where  the  displacement  increments  during  the  first  loading  cycle  (Au^  ,  Au.^) 
can  be  obtained  as  in  Fig.  4.4,  with  a  zero  initial  displacement 


u  =  u  =0 
n-1  o 


(4.16) 


Eqs.  (4.10)  and  (4.15)  give  an  upper  bound  to  n  in  the  form 

c 


n  <  n  . 
c  c+ 


iUll"Au12‘ 


(4.17) 


To  determine  a  lower  (safe)  bound  to  nc  ,  suppose  that  a  displacement 
um_^  can  be  chosen,  such  that  performing  the  operations  in  Fig.  4.4, 


Vl  +  4V  5  G 


U  +  Au,  .  . V,  >  7T 
m  (m+l)l  G 


(4.18) 


In  other  words ,  if  um_^  were  actually  the  displacement  at  the  beginning  of 
the  m**1  cycle,  m  =  nc  .  Then,  by  the  same  argument  leading  to  (4.17), 


nc  >  nc-  “  G(Au  ,-Au  .) 

ml  fim2 


(4.19) 


The  bounds  (4.17)  and  (4.19)  are  diagrammatically  indicated  in  Fig.  4.5a. 

A  value  for  um_^  can  be  obtained  in  the  following  way  (Fig.  4.6a). 
Choose  a  time  t^  »  draw  a  line  r  such  that  A  .  =  A^  »  then  perform  the 
integrations  in  Fig.  4.4  starting  from  tm1  and  the  corresponding  u^  , 
and  repeat  until 


“'W  =  Vi  +  %i  =  -• 


(4.20) 


If  (4.13)  holds. 


-(k  +  Gu)  =  -2k  <  £  at  t  =  t 


(4.21) 
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and  (4.18)  are  certainly  verified.  But  note  that  if  nc  is  large  as  assumed, 
the  inclination  of  r  must  be  very  small,  and 


m2 


=  t  .  +  T 
ml 


(4.22) 


therefore,  it  seems  fair  to  assume  that  (4.21)  holds  also  if  (4.13)  does  not, 
with  the  only  exception  of  loading  programs  of  the  type  shown  in  Fig.  4.6b. 

Even  in  this  very  particular  case,  however,  the  procedure  of  Fig.  4.6a  yields 
a  very  small  Au^  >  so  that  it  can  be  assumed  that  (4.18)  are  again  verified. 

Consider  now  the  case  in  which,  if  the  effect  of  the  axial  loads  were 
neglected,  the  conditions  of  Section  1.2  (or  the  so-called  Case  II)  would  be 
realized.  Comparing  Eqs.  (4.15)  and  (1.28),  it  can  be  noted  that  the  two  oppo¬ 
site  effects  interact;  however,  in  the  long  range,  the  instabilizing  effect 
does  prevail. 

If  the  final  residual  displacement  (neglecting  axial  loads)  un  is  much 
smaller  than  the  critical  displacement  uc  ,  it  could  be  (unsafely)  neglected 
altogether.  In  this  case,  Au^  and  Au^  introduced  in  (4.17)  must 

be  computed  starting  from  uq  =  0  and  t*  determined  as  in  Fig.  1.4. 

Conversely,  if  u  can  be  at  least  approximately  evaluated,  it  is  safe 
to  consider  it  acting  from  the  very  beginning  of  the  loading:  then,  instead 
of  nc  ,  a  smaller  number  n*  is  defined  (Fig.  4.5b)  and  bounded  between 


n* 

c- 


k 

G  ~  UR 

AulYAul2 


G  -  UR 
Auml-Aum2 


(4.23) 
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Here,  Au*^  and  Au*2  must  be  determined  with 


u  .  =  u  =  u_ 
n-1  o  R 


(4.24) 


Even  more  simply,  let  Au  be  the  increment  of  displacement  during  a  load¬ 
ing  cycle  calculated  in  the  assumption  that  the  (effective)  shear  strength  of 
the  frame  is  constantly  nil  (Fig.  4.7).  Au  is  an  upper  bound  to  the  displace¬ 
ment  of  the  frame  in  each  cycle  for  both  Cases  I  and  II,  so  that  a  safe  value 
of  nc  is  given  very  quickly  by 


n 


GAu 


(4.25) 


Noting  that 


Au  = 


't(n+m  dt 

ft  ..  f 

u  dt  =  - 

t  ,  ■ 

t  ,  • 

nl 


nl 


(n+l)l 


nl 


(t-  t  )ii  dt 
nl 


(4.26) 


is  the  first  moment  with  respect  to  t  =  tRl  of  the  two  equal  areas  A  and 
A02  4,7a  (each  taken  with  its  sign)  and  that,  because  of  (1.06),  the 

centroids  of  A^  and  Ano  are  T/2  apart,  it  is  immediate  that 


02 


a£  «  An_  J-  «  5  -- 
01  2  a  4 


(4.27) 


where  £a  is  the  average  absolute  value  of  £(t)  .  in  other  words,  the  incre¬ 
ment  of  displacement  in  a  frame  of  zero  strength  during  a  cycle  of  periodic- 
symmetric  loading,  does  not  depend  on  the  shape  of  the  loading  function,  but 
only  on  its  average  value  and  the  square  of  its  period. 

Introduction  of  (4.27)  and  (4.03)  into  (4.25)  yields 


V  =  4 


_kk. 


ZM 


g£a  T2  5  T2  ew 


(4.28) 
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Although  Eq.  (4.28)  can  only  be  expected  to  give  a  very  rough  lower  limit  on 
nc  ,  it  may  be  used  to  evaluate  the  qualitative  effect  of  some  parameters. 

Of  course ,  the  motion  of  a  rigid-plastic  frame  is  excited  only  if  the 
maximum  value  of  the  load,  ,  is  larger  than  the  frame  strength  k  .  In 

other  words 

5 

n  =  -  if  1  (4.29) 

C  K 


Thus,  any  curve  or  formula  for  n^  has  a  lower  cut-off,  that  can  be  written 
in  the  form 


t  _  r~ 


m 


(4.30) 


The  ratio  is  comprised  between  0.5  and  1  for  any  "reasonable"  loading 

function  £(t)  .  Eqs .  (4.28)  and  (4.30)  are  plotted  in  Fig.  4.8. 

Finally,  it  is  worth  noting  for  comparison  that,  after  lengthy  numerical 
calculations  of  the  responses  of  an  elastic-plastic  frame  subject  to  random 
loading  functions  and  statistical  analyses  of  the  results ,  Jennings  and  Husid 
concluded,  in  their  already  quoted  paper  [1968],  that  the  average  time  to  col¬ 
lapse  was  proportional  to  the  storey  height  and  to  the  square  of  the  ratio  be¬ 
tween  the  strength  of  the  frame  and  the  average  intensity  of  the  ground  acceler¬ 
ation;  they  did  not  investigate  the  effect  of  the  frequency  of  the  load.  Intro¬ 
ducing  the  present  symbols  and  denoting  by  C  an  appropriate  constant,  Jennings 
and  Husid' s  results  can  be  represented  by 


(n  T)  =  Ch  ( -— ) 
c  avg  X 


(4.31) 


The  qualitative  agreement  between  Eqs.  (4.28-30)  and  Eq.  (4.31)  is  remarkable: 
it  has  been  shown  [Augusti,  1969]  that  the  quantitative  agreement  is  also  fairly 


good,  and  this  would  seem  to  support  the  opportunity  for  further  researches  of 


the  type  set  forth  in  the  present  report. 
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SUMMARY  AND  CONCLUSIONS 

The  present  report  has  discussed  several  aspects  and  examples  of  behavior 
of  rigid-plastic  structures  subjected  to  dynamic  loading. 

In  particular,  in  Section  1  it  has  been  shown  that  the  response  of  strength- 
symmetric  structures  to  periodic -symmetric  loading  (in  absence  of  the  instabiliz- 
ing  effect  examined  in  Section  4)  is,  or  rapidly  tends  to,  a  periodic  motion.  In 
other  words,  the  structure  shakes  down  to  a  rigid-plastic  vibration  of  limited 
amplitude  and  the  same  period  as  the  external  agency.  The  total  accumulated 
plastic  work  (high  strain  fatigue)  becomes  the  key  factor  of  design. 

Whenever  the  response  mechanism  is  unique  (as  it  is  shown  to  be  usually  the 
case  for  plane  shear  frames  subjected  to  horizontal  loads),  very  simple  semigraph- 
ical  operations  allow  both  to  follow  the  response  cycle  by  cycle  and  to  determine 
directly  the  final  (or  steady-state)  motion. 

It  has  not  been  found  possible,  for  a  general  loading  program,  to  obtain  a 
simple  but  rational  approximation  of  the  response,  such  as  the  mode  approximation 
proposed  by  Martin  and  Symonds  [1966]  for  impulsive  loading  problems  and  recently 
extended  to  single-pulse  loads.  In  the  particular  case  of  plane  shear  frames, 
however,  even  the  multi-mechanism  response  is  governed  by  very  simple  relations. 

On  the  contrary,  the  dynamic  response  of  three-dimensional  shear  frames 
(studied  in  Section  3)  is  of  rather  cumbersome  determination  because  two  inter¬ 
acting  degrees  of  freedom  are  involved ,  except  in  some  very  special  cases .  For 
example,  the  response  to  a  periodic  load  is  given  by  a  cyclically  variable,  but 
in  most  cases  never  vanishing,  horizontal  velocity  vector.  It  is,  on  the  other 
hand,  possible  to  define  principal  directions ,  such  that  if  the  load  is  parallel 
to  one  of  them  the  response  is  one-degree-of-freedom,  and  distinguish  them  in 


stable  and  unstable. 


Finally,  in  Section  4,  the  effect  of  constant  gravity  loads  on  the  response 
of  plane  shear  frames  has  been  examined.  Firstly,  impulsively  loaded  frames 
have  been  considered,  and  a  relationship  has  been  established  between  the  total 
column  strength,  the  height  of  the  ground  storey,  and  the  momentum  that  the 
frame  can  absorb  without  collapsing. 

Then ,  it  has  been  shown  that ,  under  periodic  loads ,  plastic  deformations  of 
alternate  sign  develop,  but  the  residual  displacement  increases  with  the  loading 
cycles  until,  when  a  critical  displacement  is  reached,  the  gravity  loads  become 
able  to  cause  by  themselves  collapse  of  the  frame.  Formulae  have  been  proposed 
to  bound  the  number  of  cycles  necessary  to  reach  the  critical  displacements. 

Both  high  strain  fatigue  and  collapse  by  gravity  loads  may  be  significant  fac¬ 
tors  in  the  design  of  such  frames. 

Applications  and  numerical  examples  of  the  procedures  outlined  in  this  re¬ 
port  are  at  present  being  developed  and  will  be  included  in  later  papers. 

The  whole  treatment  presented  here  has  been  limited  to  the  consideration  of 
rigid-perfectly  plastic  behavior  and  has  therefore  the  same  limitations  of  all 
previous  works  which  have  been  based  on  the  same  assumption .  These  limitations 
will  not  be  discussed  here,  but  it  seems  appropriate  to  note  that  two  of  them 
might  be  particularly  relevant  for  the  sort  of  problems  treated  in  this  report; 
namely, 

a)  the  yield  stress  (or  yield  profile)  varies  with  the  rate  of 
deformation;  this  variation  may  not  only  cause  quantitative 
differences  (as  in  all  dynamic  plastic  problems)  but  also 
alter  the  interaction  between  the  different  degrees  of  free¬ 
dom: 

b)  the  elastic  deformations  under  a  periodic  loading  can  be  in¬ 
creased  by  resonance;  their  effect  might  be  assimilated  to 
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a  substantial  reduction  of  the  (apparent)  plastic 
strength . 

While  it  must  be  stressed  that  the  present  study  has  no  pretense  to  be 
conclusive,  and  that  the  above  problems  (and  many  others)  can  furnish  a  wealth 
of  material  for  further  investigations,  it  is  also  fair  to  remark  that  the  sim¬ 
plifications  inherent  in  the  rigid -plastic  assumption  have  allowed  to  obtain 
some  fairly  general  indications  regarding  problems  either  never  tackled  before 
or  solved  only  numerically. 
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NOTATION 


F 

f 

G 

g 

h 

I 

i 

J. 

1 

k 

L 

% 

M 

M 

o 


areas  (integrals  of  acceleration) 

axis  of  j**1  column  (in  horizontal  plan) 

axis  and  center  of  (elementary)  rotation.  Figs.  3. 1-3. 2 
centroid  and  plastic  shear  center  of  frame  in  Fig.  3.1 

total  force  (including  inertia  term) 
line  of  action  of  force  P 

parameter  measuring  the  instabilizing  effect  of  vertical 
loads  in  Eqs.  (4.01) 

acceleration  of  gravity 

height  of  storey 

intersection  of  n  and  yield  profile 
suffix  for  generic  storey,  loading  cycle,  etc. 
number  of  columns  at  i^  storey  of  frame 

strength  parameter 

total  loading  point  or  terminal  point  of  vector  £  -  u 
(Fig.  3.4  ff.) 

beam  span  (Fig.  2.1) 

bending  moment 

limit  (full  yield)  moment 


m  mass 

N  number  of  storeys  in  frame 

n  number  of  loading  cycles  (also,  suffix) 

nc  "critical"  number  of  cycles 


0 


origin  of  reference  axes  or  vectors 
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P 


p.  q 

T 


t 


t 

c 


U 

u 

u 


u 

u 


applied  force 

suffixes  denoting  particular  storeys  in  multi -storey 
frame 

period  of  time -dependent  loads 
time 

"critical  time"  or  "time  to  collapse"  (Section  4.2) 

total  response  time  for  impulsive  loading 

terminal  point  of  vector  u  (Fig.  3.6  ff.) 

displacement 

velocity 

acceleration 

ground  acceleration 


u 


R 


u 


s 

9 


u 


V. 

X 


V  . 

Ol 


V 

W 

X 

x,  y 


Au 

n 

5  =  6(t) 


£ 


m 


residual  displacement 

"steady  state"  displacement,  velocity,  acceleration 
total  shear  force  at  i ^  storey 

th 

limit  shear  force  (shear  strength)  of  i  storey 

curve  of  points  U 
vertical  loads  on  frame 

loading  point,  or  terminal  point  of  vector  5  (Fig. 
3.4  ff.) 

coordinate  axes  (also,  suffixes  for  x-  and  y-compon- 
ents  of  vectors) 

displacement  increment 

principal  direction  of  yield  profile  $ 

load  parameter  (function  of  time) 

average  absolute  value  of  5(t) 


maximum  value  of  £(t) 


-51- 


$ 

4  =  o 


yield  profile 

equation  of  yield  profile 


Notes: 

a)  Underlined  symbols  indicate  vectors  (e.g.  £  ,  u  ,  etc.) 

b)  Many  other  suffixes,  asterisks,  etc.,  have  been  often  used 
to  qualify  several  symbols. 


-52- 


REFERENCES 


Augusti,  G.,  Discussion  of  the  paper  by  Jennings  and  Husid  quoted  below; 

Journal  of  the  Engineering  Mechanics  Div.  ,  A.S.C.E.,  1969  (to 
be  published). 

Augusti,  G.,  Martin,  J.  B.,  and  O'Keeffe,  J.  D.,  "An  Approximate  Method  of 
Analysis  for  Rigid-Plastic  Pulse-Loaded  Structures."  Submitted 
for  presentation  at  the  2nd  Canadian  Congress  of  Applied  Mechan¬ 
ics;  Waterloo,  Ontario,  May  1969. 

Ballio,  G,  "Elastic -Plastic  Dynamic  Behavior  of  a  Beam-Column  Model." 

Meccanica  (Journal  of  the  Italian  Association  of  Theoretical  and 
Applied  Mechanics  AIMETA),  Milan,  Vol.  3,  No.  3,  pp.  177-186, 

Sept.  1968. 

Frick,  A.  J.,  and  Martin,  J.  B.,  "The  Plastic  Deformation  of  a  Cantilever 
Subjected  to  Biaxial  Bending  by  an  Impulsive  Load  at  its  Tip," 
Report  No.  NSF-GK1013/15 ,  Div.  of  Engineering,  Brown  Univ.,  Provi¬ 
dence,  R.  I.,  June  1968 

Harrison,  H.  B.,  "The  Ultimate  Resistance  to  Horizontal  Shear  of  a  Simple 
Space  Frame , "  Proceedings  of  the  Institution  of  Civil  Engineers, 
London,  Vol.  39,  No.  2,  pp.  313-322,  Feb.  1968. 

Housner,  G.  W.,  "Behavior  of  Structures  During  Earthquakes,"  Journal  of  the 
Engineering  Mechanics  Div.,  A.S.C.E.,  Vol.  85,  No.  EM4,  pp.  109- 
129,  Oct.  1959. 

Housner,  G.  W.,  "The  Plastic  Failure  of  Frames  During  Earthquakes,"  Proceed¬ 
ings  of  the  Second  World  Conference  on  Earthquake  Engineering, 

Tokyo  and  Kyoto,  Vol.  2,  pp.  997-1012,  1960. 

Jennings,  P.  C.,  and  Husid,  R. ,  "Collapse  of  Yielding  Structures  During  Earth¬ 
quakes,"  Journal  of  the  Engineering  Mechanics  Div.,  A.S .C.E . ,  Vol. 
94,  No.  EM5 ,  pp.  1045^1065 ,  Oct.  1968. 

Johnson,  G.  L.,  and  Martin,  J.  B.,  "The  Permanent  Deformation  of  a  Portal 
Frame  Subjected  to  a  Transverse  Impulse,"  Report  No.  BU/DTMB/12, 

Div.  of  Engineering,  Brown  Univ.,  Providence,  R.  I.,  Aug.  1967. 

Kaliszky,  S.,  "Approximate  Solutions  for  Impulsively-Loaded  Inelastic  Struc¬ 
tures  and  Continua."  Presented  orally  at  the  12th  International 
Congress  of  Applied  Mechanics,  Stanford,  Calif.,  Aug.  1968.  (Full 
text  to  be  published). 

Lee,  L.S.S.,  and  Martin,  J.  B.,  "The  Behavior  of  a  Rigid-Plastic  Column  Sub¬ 
jected  to  a  Transverse  Dynamic  Disturbance,"  Report  No.  BU/DTMB/8, 
Div.  of  Engineering,  Brown  Univ.,  Providence,  R.  I.,  April  1966. 


-53- 


Martin,  J.  B.,  "A  Note  on  the  Uniqueness  of  Solutions  for  Dynamically  Loaded 
Rigid-Plastic  and  Rigid-Viscoplastic  Continua,"  Journal  of  Applied 
Mechanics ,  A.S.M.E.,  Vol.  33,  No.  1,  pp.  207-209,  March  1966. 

Martin,  J.  B.,  "Mode  Approximations  for  Impulsively  Loaded  Structures  in  the 

Inelastic  Range,"  to  be  published  in  the  Proceedings  of  the  Interna¬ 
tional  Conf.  on  Structure,  Solid  Mechs.  and  Engrg.  Design  in  Civil 
Engineering  Materials,  Southampton,  England,  J.  Wiley,  1969. 

Martin,  J.  B.,  and  Lee,  L.S.S.,  "Approximate  Solutions  for  Impulsively  Loaded 
Elastic-Plastic  Beams,"  Journal  of  Applied  Mechanics,  A.S.M.E.,  Vol. 
35,  No.  4,  pp.  803-809,  Dec.  1968. 

Martin,  J.  B.,  and  Symonds,  P.  S.,  "Mode  Approximations  for  Impulsively-Loaded 
Rigid-Plastic  Structures,"  Journal  of  the  Engineering  Mechanics  Div., 
A.S.C.E. ,  Vol.  92,  No.  EM5 ,  Oct.  1966. 

Oien,  M.  A.,  and  Martin,  J.  B.,  "Approximate  Methods  for  Rigid-Plastic  Struc¬ 
tures  Subjected  to  Dynamic  Loading,"  Report  No.  BU/DTMB/7,  Div.  of 
Engineering,  Providence,  R.  I.,  Aug.  1965. 

Santathadapom ,  S.,  and  Chen,  W.  F.,  "Interaction  Curves  for  Sections  Under 
Combined  Biaxial  Bending  and  Axial  Force,"  Fritz  Engineering  Lab. 
Report  No.  331.3,  Dept.  Civil  Engrg.,  Lehigh  Univ.,  Bethlehem, 

Penn.,  August  1968. 

Wittrick,  W.  H.,  "Sidesway  Collapse  Modes  for  Groups  of  Vertical  Columns," 
Internat.  Jour.  Mechanical  Sciences,  Vol.  10,  pp.  549-562,  1968. 


-54- 


LIST  AND  CAPTIONS  OF  FIGURES 


Figure  1.1  -  Simple  example  of  dynamically  loaded  structure. 

Figure  1.2  -  Periodic-symmetric  loading;  Case  I. 

Figure  1.3  -  Periodic-symmetric  loading;  Case  II. 

Figure  1.4  -  Periodic-symmetric  loading;  Case  II.  Periodic  solution. 

Figure  2.1  -  Different  multiple-mechanism  behaviors  of  a  dynamically- 

loaded  simply-supported  beam. 

Figure  2.2  -  Shear  frame. 

Figure  2.3  -  Two-mechanism  motion  of  rigid-plastic  shear  frame.  (Load 

and  acceleration  diagrams.). 

Figure  2.4  -  Two-mechanism  motion  of  rigid-plastic  shear  frame.  Steady 

state  solution. 

Figure  3.1  -  Three-dimensional  shear  frame. 

Figure  3.2  -  Moment  and  shear  yield  profiles  of  j**1  column;  normality 

rule . 


Figure  3.3 

Figure  3.4 

Figure  3.5 
Figure  3.6 
Figure  3.7 
Figure  3.8 
Figure  3.9 

Figure  3.10 


Total  shear  yield  profile  and  lines  of  action  of  forces 
causing  translational  collapse  of  a  platform  supported  on 
three  columns. 


Yield  profiles  and  principal  directions  n 

Response  to  impulse  .  General  yield  profile. 

Response  to  impulse  uq  .  Particular  example . 

Motion  with  a  "strip"  yield  profile  (k  =  00 )  . 

Response  to  periodic  loading  (rectangular  yield  profile; 
both  load  components  in  Case  II). 

Steady -“State  velocities  for  loading  in  Figure  3.9. 


Figure  3.11  -  Response  to  periodic  loading  (rectangular  yield  profile; 
x-component  in  Case  I,  y-component  in  Case  II). 


Figure  3.12  -  Steady-state  velocities  for  loading  in  Figure  3.11. 
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Figure  3.13  -  Determination  of  sign  of  steady-state  movement  of  L 

and  U  . 

Figure  4.1  -  Shear  frame  under  horizontal  and  vertical  loads. 

Figure  4.2  -  Motion  of  frame  subjected  to  horizontal  loads  5(t)  and 

constant  vertical  forces. 

Figure  4.3  -  Response  of  frame  subjected  to  constant  gravity  forces 

and  a  horizontal  impulse :  a)  velocity  of  deformation 
vs.  time  t  and  initial  velocity  u  ;  b)  initial 
velocity  u^  vs.  total  response  time  t^  . 

Figure  4.4  -  Frame  subjected  to  periodic  horizontal  load  and  constant 

vertical  forces;  determination  of  n^  cycle  of  response. 

Figure  4.5  -  Bounds  to  the  critical  number  of  cycles  n  :  a)  Case  I 

b)  Case  II.  c 

Figure  4.6  -  Final  cycle  of  loading. 

Figure  4.7  -  Cycle  of  loading  with  zero  effective  strength. 

Figure  4.8  -  Approximate  lower  bound  to  the  critical  number  of  cycles. 
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